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ABSTRACT 

An understanding of stellar acoustic oscillations is attempted, using the local wave concept 
in semi-analytical calculations. The local homogeneity approximation allows to obtain simplified 
equations that can nevertheless describe the wave behavior down to the central region, as the 
gravitational potential perturbation is not neglected. Acoustic modes are calculated as classical 
standing waves in a cavity, by determining the cavity limits and the wave phases at these limits. 
The internal boundary condition is determined by a fitting with an Airy function. The external 
boundary condition is defined as the limit where the spatial variation of the background is 
important compared to the wavelength. This overall procedure is in accordance with the JWKB 
approximation. When comparing the results with numerical calculations, some drawbacks of the 
isothermal atmosphere approximation are revealed. When comparing with seismic observations 
of the Sun, possible improvements at the surface of solar models are suggested. The present 
semi-analytical method can potentially predict eigenfrequencies at the precision of ±3 /iHz in the 
range 800-5600 fiHz, for the degrees I = — 10. A numerical calculation using the same type 
of external boundary conditions could reach a global agreement with observations better than 1 
/iHz. This approach could contribute to better determine the absolute values of eigenfrequencies 
for asteroseismology. 

Subject headings: Sun: helioseismology - Sun: oscillations - Sun: atmosphere - Stars: oscillations - 
Stars: atmospheres 



1. Introduction 

Asteroseismology describes mainly the propa- 
gation of matter waves and their eigenmodes in 
stars (or planets) that are considered as cavi- 
ties. By an inverse procedure one can then de- 
duce the physical properties of the layers crossed 
by the waves. It is a very effective tool, and of- 
ten the unique one to study the internal struc- 
ture of stars. Since several decades, many in- 
tensive works were undertaken, especially for the 
Sun, to calculate precisely the oscillation modes 
and to understand the physical mechanisms in- 
volved. Different reviews are dedicated to this 
subject (Cox 1980; Unno et al. 1989; Christensen- 
Dalsgaard & Berthomieu 1991; Turck-Chieze et al. 
1993; Christensen-Dalsgaard 1998). 



Semi-analytical approaches have been devel- 
oped, essentially based on various asymptotic ap- 
proximations. Lamb (1909, 1932), Cowling (1941), 
Tassoul (1980), Leibacher & Stein (1981), Duvall 
(1982), and Deubner & Gough (1984), have ex- 
plored the first-order asymptotic approximation, 
and have helped to interpret the high-order modes 
measured for the Sun. Then with the availabil- 
ity of low-order modes, Vorontsov (1991), Lopes 
& Turck-Chieze (1994), Roxburgh & Vorontsov 
(2000), have taken into account the effect of grav- 
ity in second-order theories, allowing to better de- 
scribe the propagation in central regions of stars. 
On another hand, numerical calculations are re- 
quired to reach the high precision in frequency de- 
termination, needed for inversion of stellar internal 
layer physics. Several high accuracy oscillation 



codes already exist, offering the computation of 
eigenfrequencies and eigenfunctions (Christensen- 
Dalsgaard 1982; Unno et al. 1989). They have 
been widely used to interpret the large amount of 
solar modes observed from Earth or space. 

Just at the beginning of the spatial-asteroseismology 
era, it is worth assessing the progress of these the- 
oretical efforts in the prediction of observations. It 
is interesting to note that the agreement between 
numerical results and observations is remarkable 
in the range 1000 - 2400 /uHz where v num ~ u a b s 
is contained in ±1 /iHz. At larger frequencies 
however, the disagreement increases gradually, up 
to 20 ^Hz at 5000 ^Hz. Semi-analytical results 
agree well with numerical ones only in that high- 
frequency range, while their dispersion at lower 
frequencies stays important even when they have 
been reduced by the second order approximation. 
Moreover, the acoustic cut-off frequency is theo- 
retically predicted around 5000 fiHz, while from 
observations it is rather around 5600 fiRz (Fos- 
sat et al. 1992; Garcia et al. 1998). Those dif- 
ficulties can be overcome in the case of the Sun 
thanks to the great number of cigenmodes acces- 
sible by observations. But such difference could 
be more problematic in asteroseismology as only 
low-degree modes (I < 3) are available and the 
stellar radius and mass are less known. 

The discrepancy between calculation and ob- 
servation is generally attributed to an incorrect 
description of the outer layers of the current so- 
lar models at the surface. Several studies have 
tried to introduce a more relevant model of convec- 
tion via the introduction of a turbulent pressure 
(Kosovichev 1995). That reduces the discrepancy 
at high frequencies, but the excellent agreement 
at low frequencies is also deteriorated (Rosenthal 
et al. 1995, 1999). The role of the chromosphere 
was also invoked, but no real improvement is yet 
obtained (Wright & Thompson 1992; Vanlommel 
& Goossens 1999; Dzhalilov et al. 2000). 

With the regular progress of observational 
instruments, these remaining problems become 
more crucial: the theory-observation discrepan- 
cies at high frequencies and the inaccuracy of 
semi-analytical methods at low frequencies are up 
to 100 times larger than observational errors. So, 
we suggest to search the reasons not only in the 
solar model at the surface, but also in the eigen- 
mode calculation itself, especially about the wave 



description in the solar interior and the bound- 
ary conditions at the surface. Concepts closer to 
classical wave physics remain to be further ex- 
plored. It can be considered as a continuation of 
the works of Leibacher & Stein (1981) and Duvall 
(1982). Some first acoustic-frequency calculations 
was attempted in Nghiem (2000), and the effect of 
gravitational potential was taken into account in 
Nghiem (2003a). In this paper, we aim to further 
formalize the concepts, by putting a special accent 
on the internal and external boundary conditions. 
This work is organized as follows. In §2, a new 
p-mode calculation is developed. Its approxima- 
tions and its boundary conditions are proved to 
be in accordance with the JWKB method. In §3, 
the obtained results are first compared to present 
numerical calculations, and some drawbacks of the 
isothermal atmosphere approximation are pointed 
out. Then in §4 comparisons with observations 
suggest a way to improve solar surface descrip- 
tions, and at the same time to evaluate the poten- 
tial of the proposed method. Concluding remarks 
are presented in §5. 

In this paper, the solar model used is the seis- 
mic model (Turck-Chieze et al. 2001; Couvidat 
et al. 2003), computed with the CESAM code 
(Morel 1997). The numerical results for solar 
oscillations are obtained with the ADIPLS code 
(Christensen-Dalsgaard 1997). Acoustic mode ob- 
servations of the Sun come from the SoHO space- 
craft (MDI data from Rhodes et al. (1997); GOLF 
data from Bertello et al. (2000), and from Garcia 
et al. (2001). The graphic representation for modes 
of degrees I = to 10 follows the symbols indicated 
in Figure 1. 
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Fig. 1. — Symbols used for the graphic represen- 
tation of the I = to 10 modes in all this paper. 
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2. The Present Approach 

The present study starts from the standpoint 
that observed solar modes, due to their regular- 
ity and their precision, are related to waves in its 
most classical sense, excluding shock waves, wave 
buckets, etc. Basically eigenfrequencies are the 
signature of waves trapped in a cavity. A wave is 
the propagation of an oscillation. An oscillation 
is a continuous exchange between kinetic and po- 
tential energy around an equilibrium state. Here 
the notion of equilibrium state is the central one. 
Where this latter ends, there is no possible oscilla- 
tion, thus no more propagation, the wave becomes 
evanescent. At this point, the wave is reflected. 
That determines the turning point, i.e. the cavity 
limit. 

Inside a huge gazeous sphere maintained by 
gravity like a star, two kinds of equilibrium can 
be considered: the equilibrium between stratums 
of different pressure, and the pressure equilibrium 
inside each stratum, that is the local homogeneity. 
Every pressure perturbation generates two oscilla- 
tions: a gravity one around an equilibrium stra- 
tum, and an acoustic one around a point of local 
equilibrium. The two induced waves do not in- 
teract between them, unless their frequencies are 
close to each other. For the Sun, excepted around 
500 /zHz where coupling has to be considered, 
acoustic and gravity waves can be studied sepa- 
rately. The present approach is also applicable for 
all stars where acoustic and gravity frequencies are 
well separated. 

Let us now concentrate only on pure acoustic 
waves. By definition, such a wave can only prop- 
agate in a locally homogeneous environment, or 
more precisely homogeneous as regard to its wave- 
length. If the pressure scale height can be seen as 
representative of the homogeneity, the propaga- 
tion condition can be expressed as 

A = — < aH p{k) (1) 

where A is the wavelength, k the wavenumber, 
-ffp(k) the pressure scale height projected on the 
propagation direction k, and a is a constant to 
be determined. To be exhaustive, this coefficient 
could be considered as the first one of a Taylor 
expansion in A, H p , H p , or similar parameters. 
These coefficients should be determined once, by 
experiments on Earth, or spatial experiments, or 



theoretical studies, but they are not specific to a 
star. When a condition of the type of equation 
(1) is no longer satisfied, the wave can no longer 
propagate and is reflected. 

The reflection conditions are different, follow- 
ing the current pressure variation. If that happens 
where the pressure is increasing, for the wave, ev- 
erything is like it encounters a wall, the reflection 
gets the properties of the 'pressure wall' type. On 
the opposite case, that are the properties of the 
'pressure vacuum' type. For this latter, in the ab- 
sence of any other force, the force balance at the 
boundary requires that the perturbed pressure at 
the discontinuity is 

p' = 0. (2) 

One must have a pressure node, corresponding to 
a displacement antinode, because, as seen below, 
there is a phase lag of ir/2 between them. We will 
see that it is the case at the upper atmosphere 
limit: the last layer freely vibrates under the ac- 
tion of acoustic waves. The situation is exactly 
the same at the free end of a spring laying hori- 
zontally on a frictionless surface on Earth. Now, in 
the presence of the gravity force, the force balance 
requires that 

p' + Pa jg'dr = (3) 

where g' is the perturbed gravity, and po is the 
equilibrium density that can be brought out the 
integral thanks to the local homogeneity approx- 
imation. Note that the sign used in this context 
is opposite to that of the continuity equation. As 
seen later, this will be the situation for the radial 
mode at the internal boundary. This looks like vi- 
brations at the free end of a spring vertically hung 
in the Earth gravity field. 

From this viewpoint, the reflexion occurs at a 
given radius, and the continuity condition must be 
regarded at the two sides of this fixed limit, thus 
the perturbed pressure considered should be eule- 
rian. This is different from the perspective where 
the limit is not fixed in space but follows the dis- 
placement generated by the perturbation, where 
the lagrangian perturbed pressure is considered. 

The wave propagation and reflection are so far 
completely defined in their principles. We arc in 
the approximation context of local homogeneity, 
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or pure acoustic wave, or of local wave, i.e. the 
wave is only denned locally. It is important to 
stress that in this approximation, the environment 
is not homogeneous, but slowly variable. This 
choice is only guided by the aim of describing 
what is called a 'local-acoustic' wave. The preci- 
sion reached with this approximation is not known 
a priori, so it will be estimated on the view of its 
expressions and its numerical results, compared to 
other methods. 

2.1. Wave Propagation 

After a perturbation, the wave propagation is 
fully described by the four eulerian functions £, 
p', p', <&', which are respectively the displace- 
ment, along with the perturbations of pressure, 
density and gravitational potential. In the hy- 
potheses of adiabacity, sphericity and linearity, 
and when the time, angular, and radial parts have 
been separated, these unknown functions are gov- 
erned by the following differential equations where 
only remains the radial part (see e.g. Christensen- 
Dalsgaard (1998), p.54): 
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with r the distance to the stellar center, and I 
the mode degree. The other symbols stand for 
medium equilibrium quantities: po, po, Co, go are 
the pressure, the density, the sound speed, and 
the gravity; H p is the pressure scale height, N the 
buoyancy frequency, and Ti the adiabatic expo- 
nent. That is in fact a system of 3 differential 
equations of fourth order. Equation (6) in p' is 
added only for the sake of exhaustiveness. 

That system of differential equations can only 
be solved numerically. An analytical solution 
needs further approximations. Due to that, re- 
sults obtained by numerical codes are more pre- 
cise and realistic than those of analytical meth- 
ods. Nevertheless, the final result crucially de- 



pends on boundary conditions that are not con- 
tained in those equations. The choice of bound- 
ary conditions comes from physical considerations, 
and analytical calculations can be helpful for that, 
especially when their results are precise enough to 
be not too far from numerical ones, or in other 
words, and that is an evidence, when the approx- 
imations used are not too crude. 

In order to find an analytical solution, we pro- 
pose to use the above suggested approximation 
of local homogeneity (but not of homogeneity), 
where the derivatives of the equilibrium quantities 
are only negligible when compared to equilibrium 
quantities themselves. Only the relative deriva- 
tives, namely 1/H p , 1/H p , and thus TV 2 , can be 
neglected. This also mathematically justify a va- 
lidity condition of the type of equation (1). In that 
context, equations (4-7) become: 
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Using equations (8-11), the derivative of equation 
(9) gives 

,2 
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When searching for a spherical solution p' 
P' must verify 



p' = 0. 

(12) 
P'/r, 



d 2 P' 
dr 2 



r w 2 +47rG i o 1(1 + 1) 



P' = 0. (13) 



This typical wave equation with slowly variable co- 
efficients has a well known solution in the JWKB 
approach (see appendix A), so that 

p'(r) = ^ cos ( - [ k r dr + ip) . (14) 



By searching a solution of the same type for <£', 
equation (11) gives the particular solution 

4ttG A 
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k r dr + ip\ (15) 
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and finally equation (9) gives 



£r(r) 



rpac^k 2 



sin 



/ k r dr + tp J . 



(16) 



Thus £ r differs by n/2 in phase with p' , p' and <&'. 
Remark also that tq and r should be in a domain 
where k r > 0. Otherwise exponential solutions 
should be considered. To be exhaustive, a general 
solution has to be added to the particular solution 
but it will play a role only for radial modes 
very near the center, so we will invoke it only in 
section 2.7. 

Equations (14-16) are the expressions of a 
spherical wave that is well defined only locally. 
The amplitude A and the phase tp are not con- 
stant but almost constant. At a given radius r, 
the radial wavenumber k r is given by 



k 



2 _ u 2 +4irGp 1(1 + 1) 



(17) 



and, for a spherical wave, its horizontal component 
is 

,2 Kl + 1) 



kh — 



Thus the wave number is 



k 2 = 



u> 2 + AirGpa 



(18) 



(19) 



The amplitude of the radial displacement is fur- 
ther evaluated in appendix B. The sound speed, 
or the phase velocity of the sound wave is (in a 
propagative medium where k 2 > 0) 



co 



v k v /l + 4 7 rGp /w 2 



(20) 



which is equal to co only when 4ttGpo/u> 2 is neg- 
ligible compared to 1. 

This propagative wave solution is surpris- 
ingly simple. Discussions about its relevance are 
needed, as well in the interior than at the surface 
of the star. 

2.2. Relevance in the Stellar Interior 

In the stellar interior, a wave solution given by 
equations (14-16) are all the more simple as the 
gravitational perturbation is already taken into 
account. There is no need to develop to higher 



orders in <!>'. The latter is considered at the same 
rank than the other perturbations £, p', and p', 
like it should be naturally. The final results ob- 
tained should be more precise and more realistic. 
Indeed, as only the relative gradients of equilib- 
rium quantities arc neglected, and not the gra- 
dients themselves, the gravity that is linked to 
the pressure gradient is not neglected. Physically, 
acoustic waves have enough short wavelengths not 
to see the gravitational stratification at large scale, 
but they will progressively pass through a non con- 
stant, hydrostatic structure governed by gravity. 
The non-negligible gravitational potential and its 
perturbation are seen by the term Airgp 1 in equa- 
tion (11). 

The consequence is that the obtained accuracy 
will be the same for all ranges of I or oj (except for 
frequencies near gravity frequencies), contrary to 
restrictions of large I or u> imposed by more drastic 
asymptotic approximations. 

Another advantage is the autocoherence con- 
cerning the spherical eigenvalue 1(1 + 1) which 
will be maintained the same all along the present 
study. There is no need to replace it by (I + 0.5) 2 , 
a trick sometimes employed at low-/ (!), whose the 
artificial character is clearly pointed out in Langer 
(1937). It is also easy now to bring an answer to 
the controversial question about the fact that the 
radial mode (I = 0) reaches or not the center: due 
to its spherical nature, described by the term in 
A/r in equation (14), the wave can never touch the 
center, whatever the value of I. The non-radial 
waves reach their internal turning points where 
k r = 0, while the case of radial waves are exam- 
ined in details in section 2.7 below. 

2.3. Relevance at the stellar Surface 

At the surface, the term AnGpo can be ne- 
glected, equation (17) becomes 



1(1 + 1) 

r c 2 r 2 

and generally, for not too large I, 
kf — — 



(21) 



(22) 



which is particularly simple. The important issue 
is that k r continuously increases toward the sur- 
face. Yet it was established by current asymptotic 
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approximations that the expression of k r goes to 
zero near the surface, and that location could be 
considered as the external turning point. Does it 
mean that equation (21) is obtained in a too ap- 
proximated context, and some terms could have 
been missed? It is actually possible to derive a less 
approximated propagation equation than equation 
(13), in the following way. 

At the surface, and for not too large I, it is 
legitimate to remove the terms 2/r, 1(1 + l)/r 2 , 
and the terms in <f>' from equations (4-5), which 
become 
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Providing that 
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cldr T lPo dr + V {A0> 
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(26) 



and using equations (23-24) in the derivation of 
equation (23), one obtains: 

d 2 & i d(r lPo ) dg r 

dr 2 Tipo dr dr 

1 + ^»k = C(27> 



4 TiHpHp T lPo dr 2 
To avoid the term in d£ r /dr, introduce K as: 

1 f 1 d ^4r) (28) 



\ 2 J Txpa dr 
to obtain the propagation equation: 
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(29) 



>K = 0. 



Equation (29) is a generalization of what is derived 
in current asymptotic approximations. Indeed, 
when g is constant, the two terms after w 2 / c o can- 
cel each other out, and if further more Ti is con- 
stant, it looks more familiar (Deubncr & Gough 



1984): 



d 2 K 



dr 2 



4H 2 



2 dHp 
dr 



K = 0. (30) 



At the surface, these last equations are indis- 
putably less approximated than those obtained in 
the local homogeneity approximation. The de- 
duced k r goes to zero near the surface. But for 
all that, are those last expressions more relevant? 
Can them be used to determine the external turn- 
ing points? 

To answer, a detour by the JWKB approach is 
necessary. Following appendix A, equations (29, 
30) admit a wave-like solution only when the con- 
ditions of equation (A7) arc verified, i.e. when 
properties of the medium vary slowly. The JWKB 
approximation also predicts that in these condi- 
tions, there is no wave reflection, or at least the re- 
flections are negligible. In the opposite case, there 
is no wave-like solution and reflections will occur. 

An important consequence is that, generally, 
in order to be able to apply JWKB, the expres- 
sion of k r (r) can only contain medium variables 
at the same derivation order. For example, if the 
equilibrium pressure po appears in k r (r), then the 
derivatives dp / dr, d 2 po / dr 2 must be negligible, if 
not the conditions of equation (A7) cannot be ver- 
ified (except for very specific cases). If it con- 
tains dpo/dr, then it should contain neither p nor 
cPpo/dr 2 , d 3 p /dr 3 . 

Another consequence is that when the JWKB 
approximation is applicable, and only in this case, 
the sign of k 2 indicates the propagation or evanes- 
cent character. It is also clear that when k r (r) = 0, 
the conditions of equation (A7) are of course vio- 
lated, and there is reflection. But we should not 
forget that reflection immediately occurs when- 
ever those conditions are no more valid, even when 
k r (r) is very far from zero. And we should keep 
in mind that those reflection conditions are not 
contained in the expression of k r alone. 

It appears therefore that the k r deduced from 
equations (29) or (30) contains pa,po as well as 
their derivatives, thus is generally rapidly variable, 
thus the JWKB approach cannot be applied, and 
its zeros have no special physical sense. It just 
means that when the relative derivatives of equi- 
librium quantities cannot be neglected, we arc in 
the evanescent part, outside the resonant cavity. 
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This is coherent with the local-homogeneity anal- 
ysis. (In fact, equation (30) indicates that the 
JWKB method can nevertheless be applied in the 
very special case where both Co and H p are slowly 
variable, so that H~ 2 dH p /dr can be neglected. It 
is the case for an isothermal atmosphere which will 
be discussed further in section 3.3). 

On the contrary, in the context of local homo- 
geneity, i.e. slowly variable medium, the k r ob- 
tained in equation (17) automatically satisfies the 
JWKB application conditions. In the present case 
of spherical modes, the internal cavity limits are 
generally situated towards the center where the 
medium is really slowly variable, and the turning 
points will be determined by k r (r) = 0, or very 
close to that. Concerning the external limit, the 
cause of rapidly variable medium should rather be 
invoked. In other words, an acoustic wave turns 
back in the interior because of its own spherical 
nature, and at the surface because of an abrupt 
medium change. This internal/external dissym- 
metry can be seen on the wave trajectories (Fig. 
2). 

From that long but necessary discussion, it can 
be concluded that conceptually, the local homo- 
geneity approximation and its cavity limit crite- 
ria are nothing but a physical translation of the 
JWKB mathematical expression. 




Fig. 2. — Wave trajectories for the mode I = 
5, n = 5. Internal and external turning points are 
indicated by the arrows, which show their different 
behaviours. 



2.4. Wave Trapping 

When the acoustic wave is reflected at the ex- 
ternal and internal limits of a star, we are in the 
presence of standing waves inside a cavity, where 
only determined frequencies will be favoured. The 
sum of forward and backward waves forms stand- 
ing waves with nodes and antinodes separated by 
7r in phase 



y(r) 



A 

— cos I tot 

r 



k r dr + ip 



+ 



A 



■ COS Ldt 



ro 
r 



k r dr + ip 



= — cosiujt H — ) cos 

r 2 



k r dr + 



(31) 

where y stands for £ r , or p' ', or p 1 or $' . Fol- 
lowing Figure 3, it is clear that the phase advance 
between the turning points n and r 2 should verify 



/' 

Jri 



k r dr = (n + a)ir 



(32) 



where n is the number of nodes along r, called 
the mode order, and a is a parameter to adjust 
following the phase conditions at r± and r 2 . For 
example, if there are two antinodes at the two lim- 
its, a — 0, if there are one node and one antinodc, 
a = —0.5. If there are two nodes, a = —1. The 
convention adopted here is to count every node of 
the displacement £ r . 

The relation (32) is largely used in different ap- 
proaches for interpreting solar acoustic eigenfre- 
quencies (Duvall 1982). Even so, attention must 
be paid that r 2 is sometimes assimilated to R, the 
star radius, while it is not the case here. This 
apparently small subtility is in fact fundamental. 
Indeed, assuming r 2 = R is a crude approxima- 
tion leading to frequency errors of tens of ^Hz. 
Any attempt to compensate it, amounts to make 
a variable change while integrating (32), therefore 
to work with a wavenumber k r that is not the real 
one. So a is not the same than that employed 
here. Moreover, the phase conditions at external 
and internal turning points are completely mixed. 
The overall consequence is that one is not working 
in the real world, but in a kind of phase space. 

Here, by working with r 2 and not R, the eigen- 
mode determination becomes similar to what is 
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Fig. 3. — Acoustic standing wave with nodes and 
antinodes separated by tt in phase. Pressure and 
displacement differs by tt/2 in phase. For the pres- 
sure, the locations of the internal node, turning 
point, fitting point, and antinode are marked by 
r n , n, r f , and r a 



done for an acoustic cavity on Earth, just like 
for any musical instrument. The cigenfunctions 
projected on the radial direction are given by the 
standing wave expression (eq. 31), and the eigen- 
frequencies are determined by the Duvall expres- 
sion (32). Thus, to obtain the full solution it only 
remains to determine the internal and external 
turning points together with the phase conditions 
at these points, i.e. namely n, r 2 and a. 

2.5. External Turning Point 

Let us begin to look at the external reflection 
conditions because they are the simplest ones, for 
the phase as well as for the location. 

Following what has been discussed in section 
2, at the surface, the acoustic wave is reflected 
when a 'pressure vacuum' is encountered. For the 
wave, everything looks like the end of the mat- 
ter. Like any matter wave encountering a 'matter 
vacuum', there is an antinode at this point for the 
displacement (think e.g. about longitudinal stand- 
ing waves along a spring with a free end). That 
determines the phase. 

For the location of reflection, we must look at 
the projection of H p onto the propagation direc- 
tion. Neglecting the surface curvature (see Figure 
4), the propagation condition, equation (1), be- 
comes 

k 2 1 2tt 

k r = k r 4 > ^H P (33) 



knowing that the phase speed (eq. 20) at the sur- 
face is Co- The constant a v , v standing for vacuum, 
remains to be determined. This task concerns typ- 
ically the boundary conditions that cannot be al- 
ready included in the initial differential equations. 
Something else than the physics of wave propaga- 
tion should be considered. In order to avoid theo- 
retical results which are not totally reliable at the 
surface, we can choose the experimental Sun's cut- 
off frequency combined with the semi-empirical at- 
mosphere of Vernazza et al. (1981). According to 
Fossat et al. (1992) and Garcia et al. (1998), the 
observed Sun's cut-off frequency is 

5600 ± 150/xiiz. (34) 

In this frequency range and at the surface, equa- 
tion (33) can be simplified to 



Co a v tl p 

Figure 5 shows that, in order that the last mode 
trapped is at 5600 /iHz, a v should be: 

a v = 12.1 ±0.3 (36) 

This value could be refined by more sophisticated 
methods. For the moment, note that a change of 
a v by +0.1, will lead to a general frequency shift 
of about -2 /jHz. 



k r ' 






/k 




H p 






/ H p /cos* 







Fig. 4. — Projection of H p on the propagation 
direction k, at the surface. The surface curvature 
of the Sun has been neglected. 

Conceptually speaking, the determination of 
the cavity external limit of a star is similar to that 
of a wind musical instrument with an open end. 
There, like for the Sun, one cavity end is not ma- 
terially delimited while the cigenfrequencies are 
precisely defined. There also, the absence of an 
elaborate theory can be overcome by the use of 
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5.4x10" i | | | | 




1.0004 1.0005 1.0006 1.0007 1.0008 1.0009 

r/R 

Fig. 5. — Determination of the coefficient a v in 
equation (35). Continuous curves: 2ir/a v H p fol- 
lowing the gaz pressure of the model C of Ver- 
nazza et al. (1981). Dashed curves: ui/c. When 
a v = 12.1, the mode at 5600 pRz is the last one 
trapped. 




0.0 0.2 0.4 0.6 0.8 1.0 

r/R 



Fig. 6. — Some examples of external turning 
points determined by the intersection of the two 
curves expressed in equation (33): k 2 /k r for n = 7 
and I = 0, 1, 2, 3 (symbols described in fig.l ), and 
2ir/12.1H p (thick continuous line). 



experimental results to determine the equivalent 
external limits (Raylcigh 1896, p. 201-202). 

So all the external reflection conditions are de- 
fined. 

In Figure 6 are given some examples of external 
turning points determined by the intersection of 
the two curves expressed in equation (33), for / = 

0, 1, 2, 3 and n = 7. In such a diagram, the 
very low frequency modes I — 0, 1 and n = 1, 
2 will have their representative curves under that 
of 27r/12.1iJ p at several locations, which means 
that they hardly exist, the Sun structure being 
too inhomogeneous from their wavelength point 
of view. 

2.6. Internal Turning Point (I ^ 0) 

In the deep interior where the structure can 
be considered as locally homogeneous, there is no 
problem of boundary conditions and thus no gen- 
uine reflection point. (Except for some very low 

1, n modes, see the section above, and for the I = 
modes, see the next section). Nevertheless, the 
spherical symmetry implies that there always ex- 
ists an innermost point, called the internal turning 
point r\ , the location where the radial component 
of the wave vector vanishes. And, when studying 
only the radial part, we need to determine this in- 
ternal turning point and the phase condition there. 
For non-radial modes, i.e. with I and fc/j ^ 0, r\ is 



given by 



k r = <^> 



1(1 + 1) ~ 1(1 + 1)' 



(37) 



When the gravitational contribution is ne- 
glected (Gpo = 0), the expression (37) is reduced 
to the well known one used in the context of cur- 
rent asymptotic approximations 

„2 , ,2 

(38) 



1(1 + 1)' 



Some turning points can already be seen on 
Figure 6 where l/k r —> oo. It is also interesting 
to represent the right or left member of equation 
(37) as a function of the radius (Figure 7 for the 
Sun), so that r\ can be immediately determined 
once the frequency and the degree of a mode are 
known. For a given u>, the larger the degree, the 
less the mode goes deep. For a given degree, the 
larger the ui, the deeper the mode goes. For the 
degrees I > 5, the curves begin to merge together 
with that where Gpo = 0, which is independent 
of I. For the smaller degrees, the difference with 
the case without gravity is appreciable, even far 
from the centre, until the middle of the solar ra- 
dius, because while po decreases with the radius, 
Cg/r 2 decreases also quickly. The curve for the 
/ = 1 modes present a further special feature: the 
gravitation term becomes much more important 
compared to Cg/r 2 , so that a gap appears in the 
radiative zone, separating the curve into 2 parts, 
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and for the concerned frequencies, the outermost 
turning point r\/R is drastically reduced to 0.137 
instead of 1 for the other degrees. 




Fig. 7. — Right or left member of equation (37) 
and (65), as a function of radius. The cases of 
I = 0, 1, 2, 3, in the presence of gravity (symbols 
as defined in Figure 1) are compared to the case 
without gravity (upper thick continuous line). 

The turning points r\/R are given in Table 1 
for different I values at typical levels of uo 2 + 
1), to be compared to the case without gravity to 
appreciate the contribution of gravity for the Sun. 



Table 1: Internal turning points r\jR in the pres- 
ence of gravity, for different I at typical levels of 
u 2 /l(l + 1), compared to the case without gravity. 



1(1+1) 


10" 5 


10" 6 


310- 7 


I = 1 


0.110 


0.134 


0.137 


I = 2 


0.168 


0.388 


0.621 


I = 5 


0.198 


0.445 


0.646 


I = 10 


0.203 


0.452 


0.650 


Without G 


0.205 


0.454 


0.652 



Globally, when the gravitational potential is 
taken into account, n is smaller, especially for 
1 = 1. Gravity pulls the low / modes toward the 
centre. 

The location of the internal turning point can 
be directly determined by a mathematical expres- 
sion, while the phase condition at this point is less 
easy to find out. A discussion about H p like at the 
surface is of no use because here its projection on 



the propagation direction vanishes. Let us come 
back to equation (13) governing P' — rp' , which 
admits the JWKB standing-wave solution 



P (r) = —= cos 

y'kr 



(-/: 



k, <lr + - n ^ ) . (39) 



This solution fails of course near r = r\ where k r 
I. A well known method to bypass this difficulty 
is to develop linearly k 2 around n 



k 2 r = ki(r - n) 
and to adopt a new variable x 

x = k\^{r — n) 



(40) 



(41) 



where k\ is the derivative of k 2 in r\. Equation 
(13) becomes 



d 2 P' 
dx 2 



+ xP' = 



which admits the general solution 

P' = CAi(-x) + DBi(-x) 



(42) 



(43) 



where C and D are the constants to be deter- 
mined, Ai and Bi the Airy functions which present 
an undulatory behaviour for positive x and an ex- 
ponential behaviour for negative x. As the Ai 
function is exponentially decreasing toward the 
center while the Bi function is increasing, one must 
choose D = 0. The next step is to fit the two solu- 
tions (39) and (43) at a fit location r / . This latter 
is usually chosen so that the corresponding Xf is 
very large, in order to use the asymptotic devel- 
opment of Ai which is sinusoidal, and to deduce 
easily the wave phase. In fact, this method intro- 
duces an important lack of precision because the 
Airy solution is only valid when (40) is valid, that 
is around x = 0. We will here adopt fitting points 
77 that are a priori not far from n, corresponding 
to small values of Xf . Consider now the standing 
solution coming from the nearest antinode to r\, 
located at r a , we have (tp — ip')/2 = in equa- 
tion (39). When matching this solution and its 
derivative to the Airy solution at the the fitting 
point: 



2A 



cos / k r dr — CAi(—Xf) 



(44) 
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2A\/k r sin 



/ 



r f 

k r dr 



Ak 3/2 dK cos 
dr 



r f 

k r dr 

Ms 



-CM\-x f )k['° (45) 
the phase of the standing wave is imposed : 

rt r i dki k\ /3 M'(- Xf ) 

/ k r dr = -arctan — rrr — : ttt ~r 

J ra 4fcj? dr k r Ai(-x f ) 



(46) 

where Ai' stands for the derivative of Ai over x. As 
a is the phase lag between r\ and a displacement 
antinode, i.e. a pressure node located in r n , it is 
given by (see fig. 3) 



an 



-r-r-f'-f 



(47) 



where the terms k r dr inside the integral signs are 
implicit. This matching has only a sense when 
equation (40) can be used. Applying this latter 
and equation (41), then calculating the first and 
third integrals of equation (47), the final expres- 
sion of a can be expressed as a function of Xf : 



an 



n 

■ — +arctan • 



1 

JX5 



1 

4^7 



Ai'(-Xf) 



Ai(-xf) 



> + -x 
' (48) 

It remains now to determine the fitting point. 
It should be chosen at the farthest limit toward r\ 
where the JWKB solution is still valid, to be also 
in the validation domain of the Airy solution. As 
seen in appendix A, the JWKB method is valid as 
long as the derivative of k r is enough small com- 
pared to k r : 

(49) 



3/2 



1 



< a 



dky 2 
— <ak r ^— 

where cL IS 0j 'small' number. In other words, the 
variation of P' must be slow enough on a wave- 
length when k r goes to zero. This condition should 
be coherent with that at the external limit which is 
of the same nature (eq. 1), except that the char- 
acteristic scale height Hf to be considered here 
does not come from the environment but from the 
function P' itself. We can define a function f 



/ = 



1 



(50) 



and should use the already determined coefficient 
a v to write the condition 



< a v Hf 



(51) 



with 



H -i__^df 
f ~ fdr 



1 dkl 
Akl dr ' 



The searched condition (49) is thus 



1 dkl 



< 



2kl dr n 



(52) 



(53) 



Using equation (40), a first fitting point can then 
be obtained 



r/i = n + 



(n/2a v 



,2/3 



1/3 



(54) 



Imposing that the ratio ki/k% is enough small, the 
condition (53) is in fact not sufficient to warranty 
the validity of the standing wave solution. When I 
and n are small, the slope k\ of the function fcj? can 
be so close to zero that k^. will also remains very 
close to zero. The wave solution loose its meaning 
because the wavelength becomes too large. Either 
k\ or k%, or their product should be enough large. 
The smallest k r is determined by the largest wave- 
length that can be contained in the cavity: 



2n 



r-i - ri 



(55) 



The smallest slope k\ should be so that k r reaches 
an 'appreciable' value , let us say corresponding to 
A r < r 2 , on a 'negligible' distance, i.e. a distance 
(r — ri) small enough compared to the wavelength. 
To be coherent with above, we use again the co- 
efficient a v and the relation (1) with the inverse 
inequality sign 



A = — > a v (r - n). 



That leads to 



k\ > a v 



(2nf 

4 



(56) 



(57) 



which, together with equation (55) allow to arrive 
to the searched condition on the product fcifc^ 



(2nY 



(r 2 - ri) 2 r| 



(58) 



Using equation (40), the second fitting point can 
be obtained 



r/2 = ri + 



a.(27r) 4 



(r 2 - ri) 2 r%kl' 



(59) 
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The final fitting point can be taken as a compro- 
mise between the conditions (53) and (58), that is 
a mean value between r/i and r/2 

,, = ,, + <^! + („, 

2(r 2 - r^Vf fc 2 

Everything is now ready for the complete cal- 
culation of the I ^ modes. 

2.7. Internal Turning Point (I = 0) 

For the radial modes 



1 = 0, k h = 0, k r = k > 0, 



(61) 



there is a priori no turning point as for I ^ 0. It is 
very tempting to conclude that these modes pass 
through the star center which is a node, since due 
to symmetrical reasons it must stay immobile. In 
this case: 



ri = 0, and . 



-0.5. 



(62) 



But this induces eigenfrequencies that are lower 
than numerical or observed results by several tens 
of /iHz. In fact, no mode can touch the center: the 
linear solution of local wave cannot be applied at 
r = 0. For a spherical wave, it is obvious that the 
center is a singular point according to its expres- 
sion (31). The wave must be stopped somewhere 
before reaching the center. As the perturbed pres- 
sure p' goes to infinity, for the wave, everything 
happens as if the surrounding pressure goes to 
zero. And everything looks like the vacuum con- 
dition at the external limit, except that here the 
decreasing 'external' pressure is directly given by 
a law in 1/r. Thus H p = r, and the condition (1) 
gives the turning point n 



h'. r 



2tt 
a v r' 



(63) 



Taking into account the numerical value of a v (Eq. 
36), this relation can be written as 



AirGpo (0.51927) : 



0. 



(64) 



When neglecting the gravity, this formula is very 
close to that using by traditional approximations 
where the term 1(1 + 1) is artificially replaced by 
(I + 0.5) 2 . Therefore, the present concept of 'vac- 
uum condition' characterised by the coefficient a v , 



could be considered as the justification of this 
usual way of doing for I = 0. 

Equation (64) can also be rewritten in a general 
form like Equation (37): 



r 2 



4nGpo 
(27r/a„) 2 ~ (27r/a„) 2 



(65) 



This expression is also represented in Figure 7 giv- 
ing the internal turning point. The curve features 
an important gap in the radiative zone, even big- 
ger than for the I = 1 modes, and for the concerned 
frequencies the largest possible n/R can only go 
up to 0.035. 

The internal turning point being determined, 
let us examine the phase at this location. When- 
ever the gravity is neglected, we are in the same 
'vacuum' type condition than at the external 
boundary (Eq. 2) : the standing wave will present 
a pressure node or a displacement antinode. Thus, 
in this case 

a = 0. (66) 

In the presence of a gravity field, the condition 
(3) must be rather applied, which can be rewritten 

as 



p'(r 1 )+p (t>'(r 1 ) = 0. 



(67) 



As r\ is very near the center, the exact solution of 
4>' must be used. Indeed, resulting from Poisson's 
equation, the gravitational potential as well as its 
perturbation must remain always finite, even at 
the center. Going back to equation (11), and with 
(f)' — F' It, the equation that has to be solved is 

HL + = A * GP ' ( 68 ) 
dr 2 r 2 Cq 

To the particular solution given by equation (15), 
must be added the homogeneous solution 



F' = Cr l+1 + ^ 



(69) 



with C, D the constants to be determined. In order 
to have a negligible gravitational potential at large 
radius, C = must be chosen. And for I = 0, D 
must be chosen so that the general solution is 



$'(r) 



4ttG A 



cos 



kdr + ip 



co sip ■ 
(70) 



12 



This keeps <fi' finite at r = 0. Equation (67) be- 
comes, for the corresponding standing waves 



- ° +1 cos 
\4irGpo 



Knowing that 



f 

Jo 



-kdr + 



2 

(71) 



(72) 



with r„ the location of the closest node to ri, and 

fT\ r-r n pT\ 



an 







(73) 



o 



where the terms kdr inside the integral signs are 
implicit, the equation determining an can be de- 
duced 



2 sin an — sin 



( a7T + l 



ri 

kdr 



(~4nGp 

As generally an -C 1, this can be simplified to: 
sin J^kdr 



■- 0. 
(74) 



an ■■ 



LU 2 /AnGpn + 2 — cos J^kdr 



(75) 



As very near the center the wave number fc can be 
considered to be roughly constant, 



L 



(76) 



'0 a v 

a further simplified expression for a can be used 
sin^ 

an = , ^ — . (77) 



4-jrGpo 



+ 2 — cos 



2tt ' 



When the gravity potential is neglected (Gpo — ► 
0), this expression is coherent with equation (66). 

3. The Whole Set of Equations and Re- 
sults 

We have now every ingredient to calculate the 
eigenfrequencies. It is worth gathering here the 
whole set of equations used for the general case in 
the presence of gravity. Given the 2 parameters 
I, n which are respectively the mode degree and 
order, a system of 2 equations with 2 unknowns uj 
and a has to be solved: 



/' 



k r dr = (n + a)ir 



(78) 



and for I ^ 



a7r = -? +arctan i 



for I = 



1 Ai'(- g/ ) 
4b/ Ai(-a:/) 



sin^ 



CH7T : 



47rGp 



2 - cos # 




with 




(81) 



27T 

k r = forl^O, k r = /or / = 0, (82) 

a v r 



and ri given by 
with 

£/ is given by 



k 2 2n 



a v = 12.1. 



x f = k{ /3 (r f -n) 



(83) 
(84) 
(85) 



where k\ is the derivative of k 2 at n, and r/ is 
given by 

Equations (78-86) are the only equations 
needed for the frequency determination. This 
shows the simplicity of the present approach. 
Three functions of the equilibrium structure are 
directly involved, co(r), po{r) and H p {r). H p {r) 
is in fact only important at the surface in the 
determination of r 2 . Po( r ) contributes via the 
gravitational potential, which plays a crucial role 
toward the center, being more important for lower 
degree modes. 

One can also notice the presence of the coeffi- 
cient a v at external as wel as internal limit condi- 
tions. Determined initially by a 'vacuum' condi- 
tion at the surface, it turns to play a key role in 
the study of the validity of the JWKB solution in 
general. 
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The above equations are obtained without ne- 
glecting the gravity effect. The induced results in 
the following sections will thus present the same 
quality at low as well as at large degree or or- 
der, except when deliberately ignoring the gravity 
term. 

In the following, the results obtained with the 
above equations are exposed and commented, for 
the modes I = to 10, and frequencies in the range 
300 - 5600 /iHz, firstly for the case without gravity, 
and secondly in the general case taken into account 
the gravity. 

3.1. Results in the Case Without Gravity 

Although neglecting the gravity is a crude ap- 
proximation, it is very instructive because the cor- 
responding physical conditions look like those of 
acoustic waves on Earth, which appear more fa- 
miliar to us. Furthermore, that will allow direct 
comparisons with the existing semi-analytical ap- 
proaches where the gravity effect is neglected or 
treated at the second order. 

We employ the equations (78-86) with Gpo — ► 
0. The results are displayed in Figure 8 for the 
turning points n , r 2 , the phase shift a and the dif- 
ference v — v num between the eigenfrequencies re- 
sulting from the present work and those computed 
numerically (The numerical calculation automati- 
cally takes into account perturbation on gravity). 

The phase shift a lies in the range [-0.25, 0]. All 
the radial modes have the same a — 0, expressing 
the presence of a 'vacuum pressure'. The non- 
radial modes have a common a value only at large 
frequencies. Indeed, for these latter, when k\, the 
slope of kf ar n, is enough large, the last part of 
equation (86) is ngligible compared to the second 
one, and equation (79) gives a constant value for 
a 

a = -0.0334 (87) 

which is independent of I and n. 

The curves representing the inner turning 
points n can be compared to those presented 
in Lopes & Turck-Chieze (1994), obtained from 
solar observed modes (Libbrecht et al. 1990), with 
a formalism using L = (I + 0.5) 2 . While given by 
difffcrent formulas, the r\ of radial and non-radial 
modes seem to be very similar. 

The outer turning points r 2 feature a remark- 
able behaviour: for frequencies beyond 800 /iHz, 



r 2 depends very few on the mode order I and much 
more on frequency. Indeed, equations (81) and 
(83) show that for low /, k T and r 2 are indepen- 
dent of I. That means that a modification of the 
Sun surface will modify in the same way all the 
modes of different degree. This fact will be ex- 
ploited hereafter. This set of very grouped curves 
represents the external profile of the Sun acoustic 
cavity. One can also remark that r 2 increases very 
quickly with frequency, right beyond 2500 /iHz, 
its curve reaches already the flatten part, very 
close to r 2 /i? = 1. The results from 1000 /iHz 
can be compared to the curve obtained by Dcub- 
ner & Gough (1984) (and detailed in Christensen- 
Dalsgaard (1998)), following a general asymptotic 
expression obtained after several variable changes. 
One has otherwise keep in mind that: 

1. Confusing r 2 with R leads to frequency dis- 
placements of tens of /iHz. 

2. For larger mode degrees, r 2 depends of 
course stronger on I following equation (83). 

Concerning the comparison to numerical re- 
sults for the eigenfrequencies, a first observation 
of these raw data leads to a rather disappointing 
report: the differences reach -32/+60 /iHz. But a 
more detailed analysis shows that the differences 
to be taken into account are much smaller. Actu- 
ally, these discrepancies consist mainly in a global 
undulation independent of I, signature of a differ- 
ence, solely at the surface, between the two calcu- 
lations. As justified in the next sections, a sim- 
ple correction of surface conditions to be in ac- 
cordance with the numerical code will cancel this 
undulation. Therefore only the thickness of the 
set of these curves is to be taken into account. 
Straight away then, the classical asymptotic char- 
acter clearly appears, the differences being pro- 
gressively bigger when the frequencies go lower, 
according to the importance of Gpo in the expres- 
sion of k r (cq. 81). By discarding the results at 
very low frequencies where the coupling with grav- 
ity modes can no more be neglected, in all the rest 
between 800 and 5100 /iHz, the differences are re- 
duced to no more than ±34/iHz to ±4/iHz. Al- 
though having different r\ and a, the radial modes 
present differences with numerical results which 
are in continuity with the non-radial ones. 
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Fig. 8. — Case without gravity. Turning points, phases and frequencies for the modes I = to 10, compared 
to numerical frequencies where the gravity perturbation is taken into account. The symbols are defined in 
Figure 1. 



3.2. Results in the General Gase (With 
Gravity) 

We use the set of equations (78-86) as they are. 
The results are displayed in Figure 9. 

The external turning points r 2 are the same 
than without gravity because the gravity does not 
play any role at the surface. Every related com- 
ment in the previous section remains totally valid 
here. 

As foreseen, the internal turning points n are 
especially smaller for low-degree, low-frequency 
modes, the gravity acts like if it pulls these modes 
toward the centre. So the realted a, which de- 
pends on n, is smaller in absolute value. These 
features are particularly noticeable for the I = 
and I = 1 modes. The common value of a at large 
frequency, a=-0.0334, remains the same. The a 
of the radial modes arc now completely different, 
and are the only ones that are positive. 



These different r\ and a, together with different 
k r in the inner parts will induce different cigenfrc- 
qucncies when gravity is taken into account. One 
can moreover predict that all the frequencies will 
be lowered. These differences are larger for lower 
degree or higher frequency modes which travel 
deeper. To appreciate that numerically, some ex- 
amples of frequencies without or with gravity, are 
given in Table 2 for various degrees and orders. 
The same type of comparison for r x was given in 
Table 1. 

The curves of Figure 9 featuring the difference 
with numerical calculation show the same global 
undulation independent of I, at exactly the same 
locations in radius than with the case without 
gravity. For the same reason, only the thickness 
of this set of curves is to be considered. One can 
then conclude that between 800 and 5100 /uHz, 
the discrepancy is drastically reduced to within 
± 3 ((iHz. And the most rcmarquablc point is 
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Fig. 9. — General case in the presence of gravity. Turning points, phases and frequencies for the modes I 
to 10. The symbols are defined in Figure 1. 
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Table 2: Frequencies in /iHz without / with grav- 
ity for various degrees and orders. 





n = 2 


n= 10 


n = 30 


= 


453.2/363.7 


1570.1/1531.7 


4291.8/4273.6 


= 1 


469.6/427.4 


1629.0/1591.3 


4360.7/4340.8 


= 2 


512.2/492.7 


1680.9/1655.6 


4427.7/4410.0 


= 5 


622.1/618.8 


1810.9/1801.8 


4601.3/4590.5 


= 10 


767.3/766.4 


1998.5/1995.7 


4840.9/4835.4 



that the former asymptotic behaviour when the 
gravity was neglected, has disappeared. Even ra- 
dial and non-radial modes, which are governed by 
different physics of internal boundary condition, 
present the same profiles. This can lead to at- 
tribute a higher degree of confidence to the present 
calculations in general, and to the coefficient a v in 
particular. 



3.3. Discussion on the Outer Boundary 
Conditions 

The gravity being taken into account, it is now 
worth examining the global undulation indepen- 
dent of I (Fig. 9) characterising the frequency 
differences with numerical results. As the solar 
model used is the same, that points out a differ- 
ence in the external boundary conditions used in 
the two frequency calculations: the present semi- 
analytical calculation applies the condition of lo- 
cal inhomogeneity, while the numerical code is run 
with the isothermal atmosphere approximation. 

Although the isothermal atmosphere approxi- 
mation leads to frequencies that are closer to ob- 
servations, several non-satisfying points remain: 

• It is not very realistics, at least when com- 
pared to a semi-empirical atmosphere of 
the type of Vernazza et al. (1981), which 
presents no place with a flat temperature 
profile. No atmosphere coming from a stel- 
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lar evolution code presents any more such a 
behaviour at the surface. 

• It is not really physics, because an infinite 
isothermal atmosphere needs to be invoked, 
whose integrated mass is infinite. 

• It requires to adjust a particular atmosphere 
to the input stellar model, thus actually 
changes the initial surface conditions of the 
latter. When discrepancies with observa- 
tions are detected, it is not possible to dis- 
tinguish the part coming from the initial sur- 
face conditions, and the part coming from its 
change by the adjustment of an isothermal 
atmosphere. 

A search for boundary conditions using directly 
the input stellar structure without any modifica- 
tion, may be desirable. Then, any difference be- 
tween two sets of frequencies could be attributed 
to differences in the stellar structures considered. 
The condition proposed here (eq. 83) is only a 
suggestion in this direction, which potentially al- 
lows to recover surface properties. To prove that, 
we can try here to find out the amosphere prop- 
erties used by numerical calculation. As the wave 
solution is adjusted to that of the analytical solu- 
tion in an infinite isothermal atmosphere adjusted 
to the initial atmosphere, this double adjustment 
amounts to modify very deeply the profiles of Ti 
and Hp while keeping unchanged their product. 
So the wave propagation is not affected, only its 
outer reflexion is concerned. Therefore we have 
to find out the H p change that cancel the fre- 
quency discrepancies. In our case, this is simple to 
do because there is a one-to-one correspondence 
between v and r 2 (upper right panel of Fig. 9) 
and between r 2 and H p (Fig. 6), thus between 
v and H p . The result is shown in Fig. 10. In 
Nghiem (2003b), it was demonstrated that this 
profile is actually the H p profile when the proce- 
dure of isothermal atmosphere adjustment is per- 
formed. This profile does not have any consistency 
with the stellar structure because it results from 
a fit with a structure having very different prop- 
erties. 

With this deduced H p profile, the final fre- 
quency discrepancy is shown in Fig. 11. It remains 
large only for low frequencies < 800 ^Hz, where 
coupling with gravity modes is no more negligi- 
ble. For frequencies between 800 and 5100 ,uHz, 



the discrepancy is on the contrary reduced as ex- 
pected to less than ±3/iHz. 




Fig. 10. — Profiles of the pressure scale height at 
the surface and zoom at the external limit. Con- 
tinuous line: initial H p coming from the stellar 
evolution code. Dashcd-line: H p expected to come 
from the double adjustment of the isothermal at- 
mosphere approximation. 
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Fig. 11. — Resulting frequency discrepancy with 
numerical results for the modes I = to 10, in the 
presence of the H p profile expected to come from 
the isothermal atmosphere approximation (Fig. 
10). The symbols are defined in Fig. 1. 



4. Comparison With Observations 

Fig 12 gives the differences between frequencies 
obtained by the present study on the solar struc- 
ture coming from the evolution code, and those 
obtained by solar observations coming from the 
SoHO spacecraft. There is a global modulation 
independent of I like with numerical results but 
different in values, going down to -40 ^Hz, and 
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up to +30 /iHz. This comes indisputably from 
problems at the surface, and can be attributed to 
the fact that the solar model used does not yet 
take into account some important physical phe- 
nomenon. In Nghiem et al. (2004) it was assumed 
that it concerns exclusively a magnetic field, which 
modifies the structure by the magnetic pressure 
and the Alfven speed. Then we aim to search for 
perturbed pressure and sound speed profiles that 
cancel the frequency discrepancies. The results 
are shown in Figure 13. That finally allows to ob- 
tain the frequency differences of only ±2/iHz in 
the whole studied frequency range (Fig. 14). In 
the same way, a magnetic variation has also been 
calculated to account for the frequency variation 
related to the solar cycle. But this study has to 
be confirmed by a more detailed one, especially 
by taking into account the density change with 
the introduction of the magnetic field. 

These results does not exclude that other usual 
phenomena (turbulent pressure, non-adiabaticity, 
etc.) can reduce the theory-observation discrep- 
ancy in frequencies. Thay have only shown the ca- 
pabilities of the proposed method to recover iden- 
tified surface properties, because no a-priori hy- 
pothesis was made on the structure. If numerical 
codes can adopt the same type of boundary condi- 
tions than here, absolute frequency values could be 
calculated at a precision well better than ±l/iHz 
in the whole frequency range. That will be very 
helpful for asteroseismology in the identification of 
the modes, when only very few low-order frequen- 
cies will be available. 




v(nHz) 

Fig. 12. — Frequency differences with observa- 
tional data for the modes I = to 10. The symbols 
are defined in Fig. 1. 
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Fig. 13. — Profiles at the surface of the pressure 
and sound speed changes, expected to come from 
a magnetic field, which cancel discrepancies with 
observational frequencies. 
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Fig. 14. — Resulting frequency discrepancies with 
observational data for the modes / = 1 to 10, after 
the changes of pressure and sound speed in Fig. 
13. The symbols are defined in Fig. 1. 

5. Conclusions-Discussions 

The local wave approximation has been fully 
and coherently applied to stellar acoustic modes. 
Conceptually, it is in total accordance with the 
JWKB approach. Pratically, it consists in clearly 
determining the cavity limits and the reflection 
conditions. The standing wave nature of the 
modes are thus firmly claimed. In this simple con- 
text, a particularly reduced number of equations 
is necessary. 

The wave propagation in the innner part is nev- 
ertheless correctly represented, as the gravity ef- 
fect is directly taken into account. No develop- 
ment into higher orders, no artificial change of the 
horizontal eigenvalue 1(1 + 1) are needed. It turns 
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out that the / = modes do not touch the centre. 
Not surprisingly, the gravity plays an important 
role for low-Z modes, being able to shift frequen- 
cies by tens of fiUz. The main result is that the 
present analytical calculation is equally valid for 
low than high frequencies and degrees, only the 
coupling with gravity modes is not considered. 

At the surface, a new boundary condition is 
proposed, directly based on the propagation condi- 
tion of an acoustic wave, which corresponds to the 
JWKB validity criteria. Although that returning 
point should perhaps be refined, some first con- 
clusions can already be drawn. A coefficient like 
a v seems to be sufficient to characterize the va- 
lidity of a JWKB solution. But the main result 
is that an alternating boundary condition leads to 
frequency shifts of the same range than present 
theory-observation discrepancies. To explain the 
latter, a search for a more appropriate boundary 
condition cannot be excluded. The method pro- 
posed here has the advantage to use the atmo- 
sphere model as it is, without any modification. 
Comparisons with observations should allow to di- 
rectly recover the real atmosphere structure. Ex- 
amples are given to find out the atmosphere used 
by numerical codes in the context of isothermal at- 
mosphere approximation, and to suggest a way to 
deduce physical properties of the solar atmosphere 
when observations are considered. 

Despite the simplicity of the proposed semi- 
analytical formalism, it seems that main physi- 
cal phenomcnons involved are taken into account. 
That is why for the Sun, a precision of ±3/jHz 
can be reached in the range 800-5600 ^iHz, and 
for I = — 10. Some calculations for other stars 
show a comparable behaviour. If now a similar 
type of external boundary conditions is used in 
numerical calculations, a potential precision well 
better than 1/iHz can be expected for the whole 
frequency range. That could represent a step for- 
ward a global agreement between the three ac- 
tors involved in hclio-, astero-seismology: semi- 
analytical calculation bringing the physical under- 
standing by predicting the good trend for the re- 
sults, numerical calculation bringing the high pre- 
cision by avoiding most of the approximations, and 
observation allowing to infer real stellar structures 
when searching to cancel theory-observation dis- 
crepancies. 
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A. The JWKB Approach 

Only a brief recall of the JWKB method is proposed here, detailed historical and theoretical considerations 
are out of this paper frame. The reader can see e.g. Elmore & Heald (1969). Let y be the unknown function 
of the variable r ruled by the propagation equation 

+ «& = O (AD 
with k r a constant. The solution is of course the wave function: 

y = Aenp(±ik r r) (A2) 
with A an arbitrary constant. Consider now the case where k r is itself a function of r, the equation 

J| + fc r 2 (r)y = (A3) 

admits no known solution in the general case, and certainly no wave solution like equation (A2). A complex 
solution of the type 

y = A(r)exp[iS(r)] (A4) 
can be tried, with A, S functions of r. By separating real and imaginary parts, one obtains: 

A = CS'- 1/2 (A5) 
2 S' + 4 S' 2 



'' 2 = *?-*V + 7W (A6) 



where S',S",S"' stands for the successsive derivatives of S. When S'" and S" are negligible compared to 
S', 

S" 2 S 1 " 

-^2 << k 2 r and — « k 2 r (A7) 
and only in that case, equation (A6) can be simplified to 

5" = ±k r (A8) 

and one has a wave-like solution: 

\fk~r 



v, = -^== exp J k r dr^j (A9) 



with C an arbitrary constant. 



B. Eigenfunction, Radial Density of Kinetic Energy, and k r 

There is sometimes a confusion between the eigenfunction and the radial density of kinetic energy, or 
related quantities... Let us first of all precise that here the eigenfunctions are of the type y(r) given by Eq. 
(31), where y stands either for p', p', 0', or £ r . 

The amplitudes A/r remain to be determined. It can be done in the context of the JWKB approximation 
(eq. A9), but that is not enough detailed because some slowly varying terms can still be contained in the 
'constant' amplitude that is infact slowly variable itself. For £ r for example, physical considerations about 
kinetic energy allow to have a more detailed determination of the amplitude. For a locally plane wave, the 
kinetic energy integrated in each wavelength (one spatial oscillation) is the same. So the kinetic energy per 
unit length is proportional to the number of wavelengths per unit length, i. e. the wavenumber: 

— oc k (Bl) 
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where s is the curviligne coordinate following the wave trajectory. Now for the motion of a locally spherical 
wave projected on the radial direction, the quantity that is conserved and proportional to k r is the radial 
kinetic energy per unit of radius, or if one prefers, its flux over any sphere of radius r: 

de k 1 { d£ r \ 2 2 uj 2 2 2 



dr 2 \ <9i / 2 
So the amplitude of ^ r is proportional to 



& oc (B3) 



That is the reason why sometimes r^fp^ £ r is called the 'eigenfunction', which is in fact proportional to the 
square root of dek /dr. Figures B show that this quantity is well proportional to \/k^. Some other times, 
that is rather r^/poco £,r which is called the 'eigenfunction', and this function has consequently a constant 
amplitude along r. 
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Fig. 15. — Scaled kV" 1 (dashed-line) superimposed on scaled 'eigenfunctions' (continuous line) coming from 
the numerical code, for three different cases of L n. 
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